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The invention relates to a method of regulating or controlling a cyclically operating 
internal combustion engine using a computation model by which the cycle or portions 
of the cycle of the internal combustion engine is, or are, divided into individual parts 
and the operating condition within each cycle part is determined using measured 
values, stored and/or applied data in order to obtain actuating variables for operating 
said internal combustion engine. 

Internal combustion engines have seen many innovations in recent years such as 
turbochargers, exhaust gas recirculation, multiple injection and/or partially/fully 
variable valve timing control systems so that there has been a considerable increase 
in the number of actuating variables available for control. The possibilities resulting 
from combining the actuating variables are generally very complex and cannot be 
sufficiently ascertained using conventional global approaches such as mean value 
models or characteristics models. 

The high demands for consumption, emissions and drivability on modern internal 
combustion engines call for control concepts that cannot be carried out without the 
current status of the engine being detected. Since many of the variables required for 
control can only be measured, if at all, using expensive sensors (meaning sensors 
that are not suited for series production), there is a compelling need for novel 
computation models. 

The computing capacities within the engine control system are strongly limited, which 
places high demands on the real-time capacity of such computation models. 

DESCRIPTION OF PRIOR ART 

If at all, current methods of calculating the operating condition of an internal 
combustion engine meet the demands placed on modern control concepts with 
unsatisfactory results. The approaches used can be divided into three groups: 
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numerical methods are based on numerical integration of the processes that 
are characteristic for the cycle over the duration thereof (e.g., four strokes = 
720° crank angle). Such type methods involve complex calculations and can 
therefore not be carried out in real time under the conditions prevailing when 
used in series. 

methods based on cylinder pressure make use of the cylinder pressure curve 
measured by a suited sensor and evaluated using appropriate thermodynamic 
methods for calculating the current operating status of the engine. However, 
the sensors that are available for performing such methods are too expensive 
to be used in series or are only suited for use on the test stand. 

further known methods rely on assumptions and/or limitations that are based 
on a certain configuration of the internal combustion engine. Such type models 
are only aimed at partial functions and cannot be generalized. 

SUMMARY OF THE INVENTION 

It is the object of the invention to develop a method by which the operating status of 
an internal combustion engine can be determined readily and quickly while still with 
sufficient accuracy so as to obtain actuating variables suited for regulating or 
controlling the internal combustion engine using electronic control units (ECU) 
available for series operation. 

The solution to this object is achieved in that the computation models for the various 
individual cycle parts are based on at least partially different assumptions and/or 
have different simplifications and that the time limits of the cycle parts are at least 
partially calculated as a function of at least one variable engine operating parameter. 
The at least one variable engine operating parameter is measured or is dictated, 
depending on the operating status of the engine, by the electronic control unit (ECU) 
for example. 

The important point of the invention is that it does not only simply reduce the intervals 
between the various computations. The limits of the cycle parts are not firmly bound 
to predetermined crank angles but are made dependent on predetermined engine 
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operating parameters. The advantage that may thus be obtained is that even map 
controlled internal combustion engines with variable valve train mechanism, variable 
injection timing and the like may be mapped in a suitable manner. Appropriate 
simplifications, which permit complete analytical mapping, may be made within the 
various cycle parts, with said simplifications however not negatively affecting the 
quality of mapping as they are accurately adjusted to this part of the working cycle. 
What matters is that, within one cycle part, the operating conditions will not 
substantially change. 

If for example a cycle part performs a portion of the intake stroke that starts with the 
intake valve opening completely and ends at a point where the intake valve is 
completely closed, one takes, as a simplification for the entire cycle part, the mean of 
the intake cross section, which facilitates modeling of the gas flow. Further, for each 
cycle part, as a simplification, the piston speed is assumed to be constant by 
approximation. The error resulting from this assumption will be retroactively 
compensated later. 

The cycle parts may be defined by the complete open condition of the intake and/or 
exhaust valves, by the combustion process, by the direction of motion of the piston, 
by the compression process and/or by the expansion process. The limits of the cycle 
parts can be determined by the position of the intake and/or exhaust valves and by 
the beginning or end of the combustion process or processes. 

The solution that can be carried out for any crank angle is calculated by portions 
starting with an initial condition defined at any transition between portions of the 
cycle, the operating status being calculated in one step at the end of a portion. The 
operating status may be determined in the same way for each of the crankshaft 
angles within this portion, though. As a result thereof, the time curve of the operating 
status may also be ascertained. 

Since the processes described by comparative processes have already been defined 
analytically, more specifically algebraically, it is possible to detect the operating 
status of each cycle part in real time. 
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There is thus provided, in a further implementation of the invention, that the operating 
status at the end of the preceding cycle part be assigned to the initial conditions of 
the next cycle part. 

The operating status is at least assigned one variable from the group comprising 
torque, mass flow, in-cylinder charge condition, energy of the exhausts and heat flow 
in the cylinders. 

Depending on the operating status to be ascertained, at least one engine operating 
parameter selected from the group comprising intake pressure, intake temperature, 
gas composition in the suction pipe, exhaust pressure, exhaust temperature, 
composition of the exhaust in the exhaust elbow, parameters of the valve train 
mechanism, combustion parameters as well as general engine operating parameters 
such as engine speed and wall temperature can be calculated. To obtain this result, it 
is not necessary to measure all of the engine operating parameters for it is also 
possible to use, in parts, results obtained from algorithms. To improve the accuracy 
of the computation process, there may be provided that at least one engine operating 
parameter be determined analytically and by measurement and that computed values 
be aligned in a well known manner, with preferably at least one engine parameter 
selected from the group comprising mass flow, in-cylinder pressure, air-fuel ratio and 
torque being determined analytically and by measurement. 

In order to simplify the computation process, there is provided that the effective 
cross-sectional area of flow of the valves be approximated by a rectangular or 
stepped curve. 

With the flexible division of the cycle, the computation process is not bound to the 
type of valve train mechanism used (fixed, partially/fully variable; number of intake 
and exhaust valves). Various combustion processes (compression or spark ignition; 
number of partial combustions) only differ by the analytical solution of the portions 
performing the combustion. Computation functions independent of the configuration 
of the internal combustion engine and is affected neither by the use of pressure 
stages (compressors, turbines, and so on) nor by devices for internal or external 
exhaust gas recirculation. 
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The method includes methodology permitting computation of conditions for which 
conventional methods require numerical integration without such an integration. The 
processes involved in charge changing and combustion are generally characterized 
by time-variant parameters (e.g., valve lift, combustion history,...). These time 
variables are approximated by simplified curves (e.g., rectangular curves), which 
permits to clearly define cycle parts. The interval boundaries are flexible although 
they are a priori known by the interval definition. The cycle parts are no longer 
dependent on the time variation of actuating variables, meaning on charge changing 
and combustion history, and can be evaluated analytically. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be described in closer detail herein after with reference to the Figs 
in which: 

Fig. 1 is a schematic illustration of an internal combustion engine for 

carrying out the method of the invention, 

Fig. 2 shows a first exemplary implementation of the method in 

accordance with the invention, 

Fig. 3 shows a second exemplary implementation of the method in 

accordance with the invention, 

Fig. 4 is a valve lift diagram. 

DETAILED DESCRIPTION OF THE PREFERRED IMPLEMENTATIONS 
Example: Charge model for variable valve train mechanism 

The following assumptions and simplifications are made: 

the intake stroke is observed; the gas condition at the exhaust is the initial 
condition (in the alternative, also with exhaust stroke) 
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calculation of the charge condition (overall mass, temperature, composition, 
pressure) depending on valve timing and the current operating point of the 
engine (speed, wall temperature) for any crank angle (i.e., its curve as well). 

effective valve cross-sections are approximated by rectangular/stepped curves 

portions with different open/closed configurations of the intake/exhaust valves 
are dealt with separately. 

each portion can be calculated in one step from operating parameters and the 
end condition of the preceding portion. 

mean value model within the portion (no integration within one portion). 
The method is based on differential equations for the enthalpy variation with time of a 



cylinder: 



dHcyi/dt = Q*waii + V cy i dpcyi/dt + I H*i 



(1) 



and after conversion: 



dPcy|/dt = 1/Vcy| (- K Pcy|dV C y|/dt + (k -1) Q* wa ii + k R I Tj m*i) 



(2) 



Derivation for simplified case: 



The following simplifications are made first: 



constant piston speed: dV cy i/dt = A Q c m 

linear expression of mass flow: m* = k T ,i (pi - p cy i) 

linear expression of heat flow: Q* W aii = k w Ac y i p cy i 



(3) 
(4) 
(5) 



Substitution yields: 



dp cyl /dt = PcyiA/ C yi (-K Ao Cm + ( k -1 ) k w Acyi - k R I Tj, k T ,i) + k R/Vcyi I Pi Tj k Tti 



(6) 



wherein 
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The solution of the simple differential equation is: 

Pcyl = [Pcyl.O " P~] (V C y|A/ C y| > o) A k~ + p- (7) 

with: p- = - (k R)/(k~ A 0 Cm) I Pi Tj k T|i (8) 

\C = - K + (K -1 ) k w Acy|)/(C m Ao) - (K R)/(Cm Ac) I Tj k T ,i (9) 

The solution for the cylinder pressure consists of two parts: 

the constant pressure (negative pressure for maintaining the mass flow) 
"polytropic" for departure from the initial condition 

The solution for the entire air mass m cy i through cylinder (2) is obtained by integration 
from equation (4) 

m cy i = J I m* dt = J I k T(i (Pi - p cy i) dt (10) 

For derivation, simplifications were used that depart from real system properties and 
need therefore to be corrected retroactively: 

constant piston speed 
linear throttle equation 



- 8 - 

The points where corrections are to be made can be defined by comparing the 
approximation solution to the numerical solutions of the corresponding simple 
differential equations. 

i) Real piston speed (for linearized throttle equation) 

Substituting in the above indicated solution of equation (7) the real piston speed Cm 
for the mean piston speed, the numerical solution for low speeds can be 
approximated quite accurately. Generally, this however calls for a speed dependent 
correction simulating the lags resulting from the time variation of the piston speed. 

ii) throttle equation (for constant piston speed) 

Depending on the linearization rule k T ,j for throttle equation (4) various pressure 
differences are obtained, which are needed for maintaining the mass flow. The 
pressures, which differ while the volume is the same, result in air mass deviations. A 
correction can be made using a conversion rule for the pressure difference calculated 
for the linearized case. 

Fig. 4 depicts an example of how the effective valve cross-section is approximated by 
a mean valve cross-section. For this purpose, the effective valve lift H is 
approximated by a rectangular lift curve H m that is equal in area. For the beginning 
and the end of the cycle part, a time ti and t 2 may respectively be defined at which 
the valve lift H of the charge changing valve amounts to 10 % of the total lift. 

The internal combustion engine 1 for carrying out the method which is schematically 
illustrated in Fig. 1 comprises at least one piston 3 that reciprocates in a cylinder 2 
and defines a combustion chamber 4 provided with at least one intake manifold 5 
and at least one exhaust manifold 6 discharging therein and therefrom respectively. 
The intake manifold 5 is controlled by an intake valve 7 and the exhaust manifold 6 
by an exhaust valve 8. A fuel injection equipment 9 directly discharges into the 
combustion chamber 4. As an alternative to, or in addition to, the fuel injection 
equipment 9 an ignition equipment may discharge into the combustion chamber 4. 
The compressor member is labeled at 10, the turbine member of an exhaust gas 
turbocharger at 11. A throttle device 13 is disposed within the suction pipe 12. 
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Downstream of the turbine 11 there is provided an exhaust cleaning device 15 in the 
exhaust leg 14. Upstream of the turbine 1 1 , an exhaust gas recirculation line 16 of an 
exhaust gas recirculation 17 is connected in branching relation to the exhaust leg 14, 
said recirculation line discharging into the suction pipe 12 downstream of the 
compressor 10 and of the throttle device 13. An exhaust recirculation valve is 
indicated at 18. 

A change in the arrangement of the optional components exhaust gas recirculation 
17, compressor 10, throttle device 13, turbine 11 and exhaust cleaning device 15 will 
not influence the computation method. 

In the suction pipe 12, pressure p Ll temperature T L and/or the composition of the 
intake gas are measured. Pressure p A> temperature T A and/or composition of the 
exhaust gas are measured in the exhaust elbow of the exhaust leg 14. Further, the 
parameters of the valve train mechanism of the intake valves 7 and of the exhaust 
valves 8 are determined, namely the control times, the effective cross sectional area 
of flow of the intake valves 7 and of the exhaust valves (as a function of the valve lift 
curve). The combustion parameters, namely the control times (injection timing, 
ignition timing) and the amount of fuel are determined. Further, general engine 
operating parameters such as engine speed n and cylinder wall temperature T w are 
ascertained. Some of these operating variables can be determined algorithmically so 
that not all of the operating variables need to be actually measured. The cylinder 
pressure p cy i needs not be measured. The operating status of the internal combustion 
engine 1 is described by the following operating parameters: torque, mass flow, in- 
cylinder charge (air mass, pressure, temperature and gas composition), energy 
content of the exhaust and wall heat flow. 

For calculating the cycle of the internal combustion engine 1 in accordance with the 
present method, said cycle is divided into cycle parts 21 through 28, 31 through 38 
that are described using simplified connections and each condition within a cycle part 
21 through 28, 31 through 38 being analytically computed from the initial condition 
and the operating parameters of the respective one of the cycle parts 21 through 28, 
31 through 38. Accordingly, the numerical integration of the entire cycle is replaced 
by the combination of integrals that have been solved portionwise first. 
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The computation models are thereby based on different assumptions and/or 
comprise different simplifications. The time limits of the cycle parts 21 through 28, 31 
through 38 are calculated as a function of at least one measured engine parameter. 
An appropriate definition of the cycle parts 21 through 28, 31 through 38 is obtained 
on the basis of the position of the intake/exhaust valves 7, 8 or the sequence of the 
partial combustions. The following possibilities are thus provided: intake valve 7 
and/or exhaust valve 8 are open or a plurality of intake/exhaust valves 7, 8 are open 
at the same time; one combustion or a plurality of superposed combustions; 
compression/expansion of the gas enclosed in the cylinder. 

Fig. 2 schematically outlines a first exemplary implementation of a cycle 20 of a four- 
stroke internal combustion engine with internal exhaust gas recirculation and one 
single combustion, said cycle being divided into several cycle parts 21 through 28. 
The cycle parts 21 through 28 are characterized by the processes of combustion B, 
expansion E, opening O of the exhaust valve 8, overlapping Ol of intake valve 7 and 
exhaust valve 8, by opening I of intake valve 7 and by compression C of the gas 
within the combustion chamber 4. The cycle 20 shown in Fig. 2 comprises 
recirculating residual gas by causing the exhaust valve 8 to open again between 
intake phase I and compression phase C. 

Fig. 3 depicts a second exemplary implementation for a cycle 30 of a four-stroke 
internal combustion engine with fixed valve train mechanism, said cycle being divided 
into several cycle parts 31 through 38. In this case, the cycle 30 comprises two partial 
combustions B-i and B 2 with the cycle part 32 between the two partial combustions B^ 
and B 2 being defined as an overlapping phase Bi, 2 between the first combustion B^ 
and the second combustion B 2 . 

The method in accordance with the invention can be used as a physical charge 
model with various configurations or combustion technologies, for example both with 
a standard valve train mechanism and a partially or fully variable valve train 
mechanism and with various combustion models. Further, models for detecting the 
gas condition in the suction pipe 12 and for detecting the gas condition in the exhaust 
leg 14 can also be used. The models mentioned can be used individually or in 
combination. 



- 11 - 



Within the scope of the invention, the gas condition can also be controlled by 
selectively varying the valve timing. 

Further, combustion and exhaust gas composition with regard to C0 2 , NO Xl particles, 
and so on can be controlled by selectively varying the amount of residual gas and/or 
the combustion parameters. 

The accuracy of the method of calculation can be considerably enhanced by aligning 
the calculated parameters with measured parameters. It thus makes sense to 
compare and match the values calculated for mass flow m cy i, cylinder pressure pcyi, 
air-fuel ratio and torque with the values measured. 

The method described permits to simply determine in real time the operating 
condition for any crank angle independent of the configuration of the internal 
combustion engine 1 . 



